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ABSTRACT 

Aims. Magnetic fields play an important role during the formation and evolution of stars. Of particular interest in stellar evolution 
is what effect they have on the transport angular momentum and mixing of chemical elements along the radial direction in radiative 
regions. Current theories suggest a dynamo loop as the mechanism responsible for maintaining the magnetic field in the radiative 
zone. This loop consists of differential rotation on one side and magnetohydrodynamic instability - the so-called Tayler instability - 
on the other. However, how this might work quantitatively is still an unsettled question, largely because we do not yet understand all 
the properties of the instability in question. In this paper we explore some properties of the Tayler instability. 

Methods. We present 3D MHD simulations of purely toroidal and mixed poloidal-toroidal magnetic field configurations to study the 
behavior of the Tayler instability. For the first time the simultaneous action of rotation and magnetic diffusion are taken into account 
and the effects of a poloidal field on the dynamic evolution of unstable toroidal magnetic fields is included. 

Results. In the absence of diffusion, fast rotation (rotation rate, flu, compared to Alfven frequency, is able to suppress the 
instability when the rotation and magnetic axes are aligned and when the radial field strength gradient p < 1.5 (where p = dlnB/dlnm 
and ru is the cylindrical radius coordinate). When diffusion is included, this system turns unstable for diffusion dominated and 
marginally diffusive dominated regions. If the magnetic and rotation axes are perpendicular to each other, fix, the stabilizing effect 
induced by the Coriolis force is scale dependent and decreases with increasing wavenumber. In toroidal fields with radial field 
gradients bigger than p > 1.5, rapid rotation does not suppress the instability but instead introduces a damping factor (i>aI'2.SI\\ to the 
growth rate, in agreement with the analytic predictions. For the mixed poloidal-toroidal fields we find an unstable axisymmetric mode, 
not predicted analytically, right at the stability threshold for the non-axisymmetric modes; it has been argued that an axisymmetric 
mode is necessary for the closure of the Tayler-Spruit dynamo loop. 

Key words. Magnetohydrodynamics: instabilities, magnetic diffusion - Stars: differential rotation, radiative interiors, mixing of 
elements 


1. Introduction 

Rotation affects the structure of stars, inducing fluid motions 
such as circulation, and can enhance the mixing of chemical ele¬ 
ments. In stably stratified radiative interiors, however, fluid mo¬ 
tions are largely restricted to isobaric surfaces, and strong turbu¬ 
lent fluid motions parallel to the isobars smooth velocity gradi¬ 
ents in shells resulting in a shellular-like rotation (Zahn 1992). 
In order to smooth the radial velocity gradients, hydrodynamic 
processes appear not to be able to transport enough angular mo¬ 
mentum between neighboring shells. Consequently, it is thought 
that some process involving magnetic fields must account for any 
angular momentum transport in the radial direction. 

Even a relatively weak magnetic held is sufficient to couple 
different parts of the star and maintain a state of nearly uniform 
rotation. Eor the interior of the Sun, for example, a held of less 
than 1 gauss would be able to transmit the torque exerted by the 
solar wind through the interior (Mestel 1953). Measurements of 
the solar core’s rotation rate (Chaplin et al. 2001) show a very 
uniform rotation. This uniform rotation may be due to a mag¬ 
netic held, but this field’s origin, configuration, and strength are 
not known. By analogy with the magnetic A stars, one might 
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speculate that a ‘fossil’ magnetic held could exist in the core 
of the Sun. Since no significant net held is seen at the surface 
(averaged over the solar cycle), the radial component of such a 
fossil would, however, have to be weak - on the order of a gauss 
or less. A held weaker than this will quickly wrap up into a 
predominantly toroidal held, under the action of the remaining 
differential rotation in the core. Eventually an instability will set 
in, limiting the growth of this toroidal held. 

Near the end of the stellar lifetime, progenitors of white 
dwarfs and supernovae go through stages where the core con¬ 
tracts and spins up while the envelope expands and spins down. 
The degree of coupling between core and envelope by a magnetic 
field in this stage will affect the rotation rates of the stellar left¬ 
overs, pulsars, and white dwarfs. It is then fair to ask whether the 
rotation rate of the remnant is determined by the initial rotation 
of the progenitor, or if a secondary process must be responsible 
(Spruit & Phinney 1998; Spruit 1998). 

Tayler (1973) showed that any purely toroidal held should be 
unstable in at least some part(s) of the star. He derived stability 
conditions that are local in the meridional plane but global in the 
azimuthal direction; indeed, it is the low azimuthal m modes that 
become unstable first. Physically, one can think of the axisym¬ 
metric mode m = 0 in the following way. A purely toroidal held 
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Bip - B^(r, 6), which must be axisymmetric in order to satisfy 
V ■ B = 0, can be thought of as a collection of discrete neighbor¬ 
ing circular flux tubes. Since all stars have a very high plasma-y6, 
magnetically-induced motions are restricted to be almost incom¬ 
pressible, V ■ V 0. Now, any flux tube whose volume is held 
constant can reduce its energy by becoming shorter and fatter, 
which in the case of a flux tube encircling the axis of symmetry 
means contracting towards the axis. Obviously there is more en¬ 
ergy to be released if the field is stronger, so that if we have two 
flux tubes (of equal volume) lying at different distances from the 
symmetry axis, if the outer tube contains a greater magnetic flux 
than the inner tube, interchanging the two will release energy. 
In a stellar radiative zone, the stable stratification largely sup¬ 
presses motion in the radial direction, v • g » 0. However, to in¬ 
terchange flux tubes some radial motion is necessary; it is impos¬ 
sible for flux tubes confined to the same spherical surface to slip 
past one another. Given that the unstable m - Q mode is axisym¬ 
metric, rotation, and thus Coriolis force, has no net effect on the 
stability conditions or growth rate, as it simply causes the annu¬ 
lus to rotate around the axis of symmetry. The non-axisymmetric 
modes are also driven by the reduction of energy by letting re¬ 
gions with strong magnetic field, and therefore a strong inwards- 
pointing curvature force, move inwards, and can be thought of 
in similar terms to the fluting instability in sunspots. 

The growth timescale of the instability is on the order of 
the time taken for an Alfven wave to travel around the star on 
a toroidal held line. This timescale is very short compared to 
the evolutionary timescale of the star, e.g. about 10 years in a 
main-sequence star with a held of 1 kG. This form of instability 
is likely to be the hrst to set in as the held strength of the toroidal 
held is increased (Spruit 1999) but it is still uncertain what hap¬ 
pens next. According to a scenario developed by Spruit (2002) 
the instability could lead to a regeneration of the poloidal held 
which will eventually be transformed into a toroidal held as a 
result of the remaining differential rotation, setting up a hydro- 
magnetic dynamo loop. The balance between wrapping-up by 
differential rotation on the one hand and the destruction of the 
toroidal held by Tayler instability on the other determines the 
strength and conhguration of the held. The rate at which angular 
momentum and chemical elements are transported through the 
stellar interior depends on the resulting held strength, conhgu¬ 
ration and the details of the instability. This process has been 
implemented as a sub-grid model in stellar evolution calcula¬ 
tions by Woosley et al. (2002); Heger et al. (2003) and Maeder 
& Meynet (2003) and has been numerically reproduced in, with 
a very ideal conhguration, simulations by Braithwaite (2006a). 
More recently the process has apparently been seen in simula¬ 
tions by Rudiger et al. (2014). ' 

In the long-term, we would like to determine the type of mag¬ 
netic held that is maintained by differential rotation in a stably- 
stratihed star, under the action of magnetic instabilities, and to 
develop a quantitative theory for the transport of angular mo¬ 
mentum and chemical elements by magnetic helds in stars. This 
will require extensive simulations which include differential ro¬ 
tation and have a spherical geometry (Boldt & Braithwaite, in 
prep.) In the meantime we aim to improve our understanding 
of the Tayler instability. Progress on this front has been made 
with local simulations by Braithwaite (2006b), and has also been 


* There have been claims that this Tayler-Spruit mechanism cannot 
work because of the lack of an axisymmetric mode; it turns out, how¬ 
ever, that this is both misleading and untrue because (a) there is indeed 
a strong axisymmetric mode, as explained below, and (b) there is no 
convincing argument that an axisymmetric mode is necessary. 


investigated in global simulations by Lander & Jones (2011). 
More recently Bonanno &. Urpin (2013a,c,b) have looked at var¬ 
ious properties of the Tayler instability, in particular the effect of 
rotation, and find that rotation is essentially unable to suppress 
the instability; one of our aims in this paper is to confirm this 
result with a different method and to ascertain the growth rate of 
the instability in the presence of fast rotation. 

Another good reason to study the Tayler instability is to gain 
a better understanding of the range of magnetostatic equilibria 
available in radiative and other non-convective stars and zones. 
Since the discovery by Babcock (1947) of a strong magnetic field 
in the A star 78 Vir, strong, large-scale fields are now known to 
exist in a subset of early-type stars as well as white dwarfs and 
neutron stars (see e.g. Braithwaite &. Spruit subm. for a re¬ 
view). These fields, known in the literature as fossil fields, are 
the relics of some previous stage of evolution, and have found 
their way into a stable equilibrium, evolving only on the very 
long diffusive timescale. An arbitrary initial magnetic held is 
known to relax in a non-convective star into such an equilibrium 
(Braithwaite & Spruit 2004; Braithwaite & Nordlund 2006), and 
there are apparently a large range of equilibria available, includ¬ 
ing both approximately axisymmetric ’twisted torus’ conhgu- 
rations as well as more complex non-axisymmetric conhgura- 
tions (Braithwaite 2008). Amongst the axisymmetric equilib¬ 
ria, it is an interesting question to ask in what ratio the poloidal 
and toroidal held components must be present to ensure stabil¬ 
ity, since both are unstable on their own. It seems that a rel¬ 
atively weak poloidal held can stabilize a toroidal held (Braith¬ 
waite 2009; Akgiin et al. 2013) but the detail of how this works is 
still uncertain. One reason why the ratio of the two is particularly 
interesting is that a star with a predominantly poloidal held is 
oblate in shape, but if the held is predominantly toroidal, the star 
is prolate. Fast-spinning highly-magnetized neutron stars (the 
so-called millisecond magnetars) with predominantly toroidal 
magnetic helds might Hip over until the magnetic and rotation 
axes are orthogonal, the damping of free precession reducing the 
kinetic energy by maximizing the moment of inertia about the 
rotation axis. Such stars would emit gravitational waves poten¬ 
tially detectable with the next generation of detectors. 

In this paper we use numerical methods to investigate the 
linear properties of the Tayler instability, much of which is still 
uncertain. We provide a useful check on the analytical results, 
which may have missed some aspects of the mechanism. We use 
similar methods to those used by Braithwaite (2006b) to inves¬ 
tigate the Tayler instability, but we investigate different aspects 
of the instability. In the next section, we look at what is known 
about the instability, before describing our numerical model in 
section 3, presenting our results in section 4 and summarizing in 
section 5. 


2. Stability of toroidal fields: analytic results 

Tayler (1973) derived the stability conditions of purely toroidal 
helds in radiative stellar interiors, ignoring the effects of vis¬ 
cosity and of thermal and, magnetic diffusion. The occurrence 
of instability for the adiabatic case is independent of the held 
strength but only depends on the held conhguration. However, 
when the instability is already present the growth rate depends 
on some measure of the held strength and the stratihcation. The 
held strength is characterized by the Alfven frequency, which 
can be thought as the inverse of the timescale taken by an Alfven 
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wave to go around the star along a toroidal field line, defined as 


Va _ 

Rt 


( 1 ) 


where Va = B^/ is the Alfven velocity in a toroidal field 
line and is the radius of the star. The stratification is charac¬ 
terized by the Brunt-Vaisala or the buoyancy frequency, N, and 
radiative zones are always in the regime 


N » a>A,<p- 


( 2 ) 


This is the same as saying that the thermal energy density is 
much bigger than the magnetic energy density, i.e. a large 
plasma- p = P/SttB^ parameter, which is true in all stars. 

Following Tayler, we make the analysis easier by using 
cylindrical coordinates {m,(p,z). We now consider an axisym- 
metric, azimuthal toroidal field independent of height z, R - 
B(r(T)e^, and look at the region around the axis. The important 
parameter is the radial field gradient p: 


P = 


dlogB 
d log nr 


as nr ^ 0. 


( 3 ) 


Now, if we perturb a toroidal field in equilibrium, the result¬ 
ing displacement can be described by 

^ gi(m0+fe)+o-r 


where the variation in the nr direction is not included as its 
length scale is much larger. Although displacements are pre¬ 
dominantly horizontal because of the vertical stratification, the 
incompressible continuity equation (whose use is justified by the 
very high plasma- P) shows that a horizontal displacement is 
necessarily accompanied by a radial displacement on the order 
of ~ ^hlrPh, where 1;, and Ir are horizontal and radial length- 
scales, respectively, so that near the axis = 1/k. The shape of 
some of the unstable modes is shown in figure 1. Note that the 
m - 0 mode is also known as the sausage mode and the m - 1 
mode as the kink mode. 

In this figure we can observe why perturbations like (4) in a 
purely toroidal field are unstable. Regions where magnetic field 
is pinched have higher magnetic pressure forcing the fluid to dis¬ 
place opposite to the pressure gradient, leading to a runaway pro¬ 
cess. Figure 1 shows the equilibrium configuration and the shape 
of the azimuthal m - 0 and m — \ unstable modes. In a stellar 
interior we often expect a vertical stratification, and vertical dis¬ 
placements experience a damping effect because they have to do 
work against the stratification. Flowever, given that the m - 1 
mode mostly drives azimuthal displacements, the damping ef¬ 
fect of the stable stratification is very weak. 

Tayler (1957) showed that for the instability to develop the 
field must satisfy 

2 

m 

p > -^— 1 {m + Q) and p > 1 (m = 0). (5) 

These instability conditions are valid only for non-viscous, non¬ 
rotating, non-diffusive, radiative stellar interiors. Examining this 
stability condition (5) we observe that for a p = 1 field, the m - 1 
mode is the only unstable mode, the m - 0 and 2 modes are 
marginally stable and the other modes are stable. For a some¬ 
what steeper radial field gradient of p = 2, modes m = 0,1,2 are 
expected to be unstable. 



Fig. 1. Physical form of the different unstable azimuthal m modes 
present in a plasma column with a purely toroidal field. Below each 
mode the cross section of the co-moving surface is shown along with 
the cross section of the equilibrium state. The only axisymmetric mode 
is m = 0 and the only mode that preserves the column’s cross section is 
the m = 1 mode. 


To study the growth of the instability we define the local 
toroidal^ Alfven frequency as 

WA,0 = VA,0/ra- =-(6) 

m Anp 

where Ma,(P is the toroidal Alfven speed and m the distance to the 
axis, we expect the growth rate of the instability to be 


O- ~ tUA,0. (7) 

Given that toroidal magnetic fields in stars are buried under the 
surface, there is no direct observational evidence of their config¬ 
uration, therefore we have no knowledge of the radial magnetic 
field gradient. Even less understood is the configuration of a 
magnetic field inside a differentially rotating star. However we 
know that towards the magnetic axis the field must tend to zero, 
and to avoid singularities the radial gradient p > 1 for ur —» 0, 
and for this reason in this paper we start by studying the p = 1 
case. Since it is expected that there is a qualitative threshold at 
p = 3/2 (see section 2.1), we then also study the case of a p = 2 
field gradient. 

The unstable growth rate in equation 7 is valid only for 
purely toroidal, non-rotating, non-stratified and non-diffusive 
cases. We now discuss the behavior of the Tayler instability 
for situations where these very ideal conditions are relaxed and 

^ For a purely toroidal field the toroidal and total (that is, replacing 
by B in the definition) Alfven frequencies are the same; however for 
mixed poloidal-toroidal fields this distinction is necessary. 
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discuss what happens when magnetic diffusion is included, rota¬ 
tion is present and the magnetic field configuration is not purely 
toroidal anymore. By doing so, the stability conditions and 
growth rates are modified, and furthermore growth rates develop 
length scale dependencies. A stable stratification has major ef¬ 
fects on the behavior of the instability, this is however outside 
the scope of this paper and is already discussed in Braithwaite 
(2006b). 

2.1. Effects of rotation 

All stars are rotating and a big fraction rotates at very high speeds 
(Collins 1963; Ramfrez-Agudelo et al. 2013). 

2.1.1. Rotation and magnetic axes aligned 

Acheson (1978); Pitts 8c Tayler (1985) and Spruit (1999) showed 
that both the stability conditions and the unstable growth rate are 
affected by fast rotation. When the rotation frequency is higher 
than the Alfven frequency (Q|| > wa,^) and the rotation and mag¬ 
netic axes are aligned, the stability condition is modified from (5) 
to 

2 

p > — - 1-1 {m + f)) and P > 1 (m-0) (8) 

so that the condition of stability against all m modes is p < 312 
rather than p < -1/2 as before, and furthermore, the growth 
rate when any unstable mode is present is damped by a factor of 
wa,(a/^^|| (see Acheson 1978; Spruit 1999). 

2.1.2. Rotation and magnetic axes not aligned 

It is not always the case that the magnetic axis, the rotation axis 
and the gravity are aligned, as in the cases looked at above. 
Figure 2 (upper panel) shows the poloidal field lines in an ax- 
isymmetric configuration, and a toroidal-field-like configuration 
is observed around the neutral field line around the star, where 
the poloidal field vanishes. This configuration suffers from the 
same kind of instability (Tayler instability when considered lo¬ 
cally), the difference being that the magnetic axis (in this case 
the neutral line) is perpendicular to gravity, but because of the 
switch in geometry, it is now the poloidal field that drives the 
instability, instead of the toroidal field as before; during local 
considerations it therefore makes sense to switch the definitions 
of the two words. Also note that in the presence of rotation, the 
magnetic and rotation axis wont be alligned somewhere in the 
star. Kink-like instabilities appear with displacements perpen¬ 
dicular to the direction of stratification, v ■ g 0, as appreciated 
in figure 2 (lower panel). Displacements are mainly in the di¬ 
rection of the arrows, but, as discussed in section 2, there are 
also necessarily smaller displacements in the direction parallel 
to the neutral line . The ratio between the displacements in these 
two directions should be equal to the ratio of length scales in the 
same two directions, as we can see from the constraint V • v 0, 
a consequence of the high plasma-/?. 

The Coriolis force affects motions perpendicular to the ro¬ 
tation axis; in the case where the magnetic and rotation axes 
are aligned, unstable displacements are mainly perpendicular 
to the rotation and magnetic axes, and stabilization in the case 
-1/2 < p < 3/2 is achieved when Qy > wa.^- When the mag¬ 
netic and rotation axes are not aligned, the amplitude of the dis¬ 
placements perpendicular to the rotation axis decrease with in¬ 
creasing inclination angle, and the extreme case occurs when the 




Fig. 2. Above: Mixed poloidal-toroidal field configuration. If the 
toroidal field is too weak or indeed absent, the poloidal field is sub¬ 
ject to an instability. Below: Left, an equilibrium configuration with 
poloidal field loops threaded by the neutral line, and on the right-hand- 
side the kink-like instability present around the neutral line. Unstable 
displacements are restricted to be approximately perpendicular to the 
stratification. Figure taken from Braithwaite 2007. 

rotation and magnetic axes are perpendicular to each other. Then 
naively we think that in order to stabilize the perpendicular rota¬ 
tor, a faster rotation rate is required, and as we assumed that the 
fluid is almost incompressible, V ■ v a; 0, we derive the stability 
condition 

> coA.^mok, (9) 

where Qj^ is the perpendicular rotation rate with respect to the 
magnetic field axis, mo is the typical length scale of the mo¬ 
tion perpendicular to the magnetic field and k is the wavenum¬ 
ber in the direction of the magnetic axis (here the neutral line) as 
defined in (4). Note that the new perpendicular rotating stabil¬ 
ity conditions now dependent on the wavenumber, unlike in the 
aligned rotator case. 

This effect can be explained if we imagine the plasma col¬ 
umn to be a series of disks subject to pressure on top and below. 
When the disks displace perpendicular to the magnetic axis they 
suffer a deformation fattening at one side and shrinking at the 
other. The magnitude of the deformation depends on the size 
of the disks, as shown in figure 3, where disk of different sizes 
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to the equilibrium energetically impossible; as magnetic energy 
is turned into heat the system moves irreversibly away from its 
original equilibrium. 

Magnetic diffusion is scale dependent, and therefore when 
included in the analysis, the unstable growth rate of the Tayler in¬ 
stability becomes scale dependent as well (see e.g. Spruit 1999). 
For a given magnetic diffusivity, rj, the diffusion rate is given by, 
rjk^, where k is the wavenumber. At first glance one expects this 
to impose an upper limit for the unstable wavenumbers that are 
able to grow; perturbations at a larger wavenumbers should be 
smoothed away faster than they can grow. We have stability if 

rjk^ > cr. (10) 

Therefore for a given wavenumber we can define a critical diffu¬ 
sivity 

^crit ~ 5 (11) 

From now on, we will express the different diffusive environ¬ 
ments as fractions of the critical magnetic diffusivity 

, ( 12 ) 

ricrit 

where one can always move back and forth between fraction of 
critical magnetic diffusivity and the wavenumber of interest with 
the following relation 


Fig. 3. Representations of the first two vertical n modes of the unstable 
;« = 1 mode. The figure on top shows the purely toroidal magnetic field 
configuration, where the orientation of the magnetic symmetry axis B 
is perpendicular to the rotation axis fl,,. It also shows what an u = 1 
and n = 2 disk looks like as the red and blue box, respectively. The 
two figures below show the slipping disks for the n = 1 (left) and n = 2 
(right) cases. The magnitude of the deformation Al±, and the defor¬ 
mation velocity Avj^, suffered by the disks is different depending on 
the unstable mode. Small n modes present large deformations, whereas 
high n modes show small deformations, A/j^(n = 1) > A/j^(n = 2) and 
Avx(n = 1) > Avx(n = 2). 

correspond to different vertical ii modes. This is simply a con¬ 
sequence of the constraint V ■ v a; 0, which can be written as 
kVj ~ The Coriolis force affects motion in the vertical 

direction; Greater deformations, in other words small n modes, 
feel a stronger stabilizing effect. Note that the Coriolis force also 
affects motions perpendicular to both the magnetic and rotation 
axis, but this motion, at least for the m = 1 mode, can be min¬ 
imized by combining left- and right-handed spiral modes with 
the right phase. 

2.2. Effects of magnetic diffusion 

Until now we have only considered adiabatic perturbations, that 
is ignoring the effects of thermal and magnetic diffusion, as well 
as viscosity. We now relax the flux freezing condition and allow 
the field lines to slip through the plasma. 

Pitts & Tayler (1985) and Spruit (1999) suggested that in the 
presence of any diffusive process, thermal or magnetic, the in¬ 
stability condition in fast rotators (8) would recover its original 
form given in condition (5) (see Spruit 1999, Appendix). The 
reason for this can be thought of in terms of energy dissipation; 
rotation tends to give stability because the Coriolis force, which 
acts perpendicularly to the velocity and has consequently no ef¬ 
fect on the energies, drives the fluid back into its equilibrium po¬ 
sition; diffusion damps these epicycle motions, making a return 


where rj is the magnetic diffusivity of the system and k - 1 /4 is 
the vertical wavenumber. 

2.3. Mixed poioidai-toroidai configurations 

The addition of a poloidal field is expected to modify or suppress 
the instability. If one imagines the toroidal and poloidal fields 
to operate independently of each other (in some vague sense) 
then some of the energy liberated from the toroidal field must 
be used to bend the poloidal field lines. If the poloidal field is 
strong, the instability might be suppressed altogether. However, 
the effect will depend on the wavenumber of the mode; whilst 
all wavenumbers release the same energy from the toroidal field, 
higher wavenumbers have to do more work against the poloidal 
field. This gives rise to a wavenumber stability threshold. 

In general, if an instability displaces a fluid element hori¬ 
zontally a distance it does so by means of a force per unit 
mass of Fh - cr^^h where cr is the growth rate. In the case of 
the Tayler instability, cr » as discussed at the beginning 
of section 2 (equation 1). Imagine, in the cylindrical geometry 
used above, that there is a poloidal field in the z direction. The 
displacement caused by the instability drags the poloidal field 
with it, which reacts against this bending with a force per unit 
mass of B^Bj^kj ^j4-np » where k is the radial 

(vertical) wavenumber. For a perturbed displacement to grow it 
is required that the instability force is bigger than the restoring 
force by bending the poloidal field line, which gives a condition 
on the unstable wavenumbers 

k<—, (14) 

Va,z 

where va.^ = B^/ ^Anp is the vertical Alfven speed. See Spruit 
(2002) and Braithwaite (2009) for a fuller discussion of this 
topic. 
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Fig. 4. Predominantly toroidal magnetic field configuration in the 
interior of a star. The box shows a zoomed-in region containing the 
toroidal field lines around the magnetic axis and bold arrows showing 
the direction of the parallel and perpendicular rotation axes. 


3. Model 

To determine the effects of rotation and magnetic diffusion on 
the Tayler instability, local simulations of a radiative stellar in¬ 
terior are performed. For our local simulations we select a box 
containing the magnetic axis of the star, located at the center of 
the box, as shown in figure 4. Global simulations of the stellar 
radiative region are not suitable at this stage because of reso¬ 
lution limitations to accurately resolve small scales. However, 
global simulations might eventually be necessary to investigate 
the non-linear behavior in a differentially rotating star, or rather, 
global in two dimensions and local in the radial direction. 

The equation of state used is the ideal equation of state, in 
units such that the molar gas constant divided by the molecular 
mass is unity 

P - pT and e - —. (15) 

7-1 


To trace the evolution of the system we solve the standard 
MHD equations, momentum equation: 


and the specific internal energy: 

^ =-V-(ev)-FV-v (19) 

ot 


where p is the density, v the velocity, B the magnetic field, 77 the 
magnetic diffusivity, e the internal energy, g the gravity, P gas 
pressure, and T gas temperature. The ratio of specific heats 7 is 
set to 5/3. We note that in the momentum equation (16) we ne¬ 
glected viscosity and the Coriolis term is explicitly included. In 
the induction equation (18) magnetic diffusion 77 is spatially con¬ 
stant. The energy equation neglects viscous and Ohmic heating 
and thermal conduction, since we are not interested at this stage 
in viscous and thermal effects. We do not include the effects of 
gravity in this study. 


3.1. Numerical code 

We use the stagger code, a parallel 3D MHD code originally 
developed by (Nordlund & Gaalsgard 1995), and continuously 
improved over the years by its user community. The Stagger- 
code uses a sixth-order explicit finite-difference scheme for 
numerical derivatives and the corresponding fifth- order inter¬ 
polation scheme. The solution of the hydrodynamic equations 
is advanced in time using an explicit third-order Runge-Kutta 
integration method (Williamson 1980). The code operates in a 
staggered, Eulerian, Cartesian mesh, which seem to be a better 
implementation than cylindrical and spherical coordinates, as 
a Cartesian mesh avoids the coordinate singularity in the axis 
present in cylindrical and spherical grid-based codes. Another 
advantage of a Cartesian based code is reflected in the simplicity 
of the equations with respect to another coordinate system. 
Cartesian coordinates also have their own disadvantages such as 
the noise introduced by fitting a cylindrical shape, e.g. torus, in 
squared boxes, or the extra computational time cost due to the 
un-used corners of the box. 

For the simulations of purely toroidal fields a resolution of 
48^ in the azimuthal direction and 32 in the vertical direction was 
used, for the mixed poloidal-toroidal fields the resolution was 
48^. Test cases with a resolution of 96^ were performed in order 
to compare the results with the low resolution simulations and 
no significant differences were observed. Periodic boundaries 
are used in all simulations; see below. 

The code also contains highly localized diffusivities that pre¬ 
vent the growth of wavelengths close to the Nyquist frequency. 
We however decided to switch off this artificial diffusive terms 
in order to avoid non-physical effects. This was possible given 
the clean initial conditions, where the development and growth 
of the Tayler instability occurs at very early times when no arti¬ 
ficial stabilization is required. 


= -V ■ (pvv) - VP -I- —V X B X B - gp H- 2pv x Q, (16) 
ot A-n 


continuity equation. 


3.2. Initial conditions 

The initial conditions correspond to the equilibrium configura¬ 
tion of a plasma containing a magnetic field of the form: 


dt 


-V ■ p\. 


(17) B = B^(ziT)e 0 -H 


( 20 ) 


induction equation, 

/9R 

— = V X [(v X B) - 77 V X B], 


The magnetic configuration in this paper has no radial compo¬ 
nent, Bjjj - 0 , and is a combination of an azimuthal, B^, and 
a vertical, B-, components with no dependence on height z or 
(18) angle </>. 
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Fig. 5. a) Initial purely toroidal magnetic configuration. The region 
of interest is between 0 < cr < cjo given that the field strength falls 
off at greater m. Colored areas represent surfaces of equal distance 
to the axis. Field lines are also plotted; their density is related to the 
increase in the field strength proportional to the distance from the axis, 
b) Mixed poloidal-toroidal field configuration. The field lines show a 
spiral morphology, for this particular case (with p = \) two spirals turns 
are completed in the box, i.e. toroidal Alfven frequency is double the 
vertical Alfven frequency 


The initial conditions for the thermodynamic variables are 
determined by the equilibrium solutions found by setting djdt - 
V = 7/ = 0 in equations (16) to (19), with the initial magnetic field 
given by equations (20, 22 and 21). No height dependence is 
considered, therefore pressure, density and temperature depend 
only on the distance from the axis. For a full derivation of these 
stability configurations see Braithwaite 2006b (section 3.3). 

Four different cases for two different field configurations are 
simulated; 

1 rotating diffusive star, 

2 perpendicular rotator, 

3 rotating star with a quadratic radial field gradient; p -2, and 

4 mixed poloidal-toroidal field. 

Case 1: Rotating diffusive star with linear radial field gradient, 
p - From (21) we see that the toroidal Alfven frequency 
varies as wa at low bt, so that it is approximately 

independent of location in the inner part of the computational 
box. The magnetic and rotation axes are aligned. We explore 
the parameter space of rotation-only, magnetic diffusion-only 
and the simultaneous action of rotation and diffusion. 

Case 2: Perpendicular rotator with linear radial field gradient, 
p - 1. Rotation and magnetic axes are perpendicular. We 
explore the rotation rate parameter space. 

Case 3: Rotating star with a quadratic radial field gradient, 
p - 2. The Alfven frequency increases linearly with the distance 
to the axis. For this case rotation and magnetic axes are aligned. 
We explore the rotation rate parameter space. 

Case 4: Mixed poloidal-toroidal configuration. The poloidal 
field is constant in the box, therefore the vertical Alfven 
frequency is also constant. The toroidal field has a linear radial 
field gradient, p = 1, thus the azimuthal Alfven frequency 
is also constant in the inner part of the box. We explore the 
poloidal-to-toroidal field strength ratio parameter space. 


The azimuthal component is an axisymmetric torus defined 
by 

B^Cct) = (21) 

\v7ol 

where btq is a constant length, set to one quarter of the horizontal 
size of the box, p is the radial field gradient, set to either 1 or 2 
in the simulations, and Bo ,0 is a constant. The exponential term 
exp(BT“/BTQ) shapes the field such that it behaves as B oc m'’ near 
the axis and vanishes farther away so as to avoid problems at the 
boundaries. The local toroidal Alfven frequency, as defined in 
( 6 ), is ua, 4 i = B^/bt ^jAnp. The maximum initial field strength 
in the box corresponds to a minimum plasma beta parameter of 
P - 10 , i.e. thermal pressure dominates over magnetic pressure. 

The initial poloidal field is implemented as a uniform field in 
the vertical direction 

B, - Bo,, (22) 

for this field configuration the vertical Alfven frequency is also 
constant everywhere, ua,z - VA,,( 2 ; 7 r/L,), where va,, = B,/ yjAnp 
is the vertical Alfven velocity and L, is the vertical extent of the 
box. 

A visual example of the field configuration, for purely 
toroidal and mixed poloidal-toroidal fields, is shown in figure 
5. 


Given that the initial configurations correspond to equilib¬ 
rium solutions, an initial perturbation is included to trigger in¬ 
stability. 


3.3. Initial perturbation 

The initial conditions in section 3.2 correspond to a perfect hy¬ 
drostatic equilibrium, in order to see any instability grow we 
seed tiny perturbations in the velocity field. Given the cylin¬ 
drical symmetry of this problem, we can perturb the different 
azimuthal, m, and vertical, ii, modes independently. The shape 
and description of the perturbation used here are the same as 
used by Braithwaite (2006b, sec. 3.5). The perturbation has its 
maximum at bt = 0 and decreases with increasing bt, almost 
vanishing at bt = btq. There is no perturbation to the vertical 
velocity component. 

The initial perturbation has the form 




m=0,n 


5 BT' 


■ 2 \ 


^0^ 




mi:0,n 


mo 


ZJ- 


Tn~c 


COS(fe)en,H- 
[cos(fe - m(p)ej;, 


+ 


sin(fe - m(p)e^]. 


The amplitude coefficients of the perturbed velocity field 
are normalized to the maximum Alfven velocity in the box: 
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V"„/vA,max = 1 X 10“*. Such a value is large enough to trig¬ 
ger the Tayler instability and small enough to let us follow the 
linear growth phase for several growth timescales. 

3.4. Boundary conditions 

The box represents a small volume of a radiative zone in a stellar 
interior. Figure 4 is a cartoon of the box location in star with 
respect to the magnetic field axis. 

The difference in length scales between the instability and 
the radius of the star are huge. The instability is many orders of 
magnitude smaller^ than the radius of the star (4 «: R,), these 
instabilities are local in the radial direction and global in the az¬ 
imuthal direction, therefore different shells of the star are in ef¬ 
fect disconnected from each other. 

Periodic boundaries are the best option for our simulation, 
given that we want to avoid the selection of special locations 
in the stellar interior and to prevent additional effects such as 
current sheets or shear. 



Fig. 6. Amplitude of the first four m modes as a function of time, ex¬ 
tracted from the magnetic field 6B(t), and normalized to the wavelength 
corresponding to the n = 1 vertical mode. The azimuthal m = I mode, 
solid red line, is the only unstable mode. 


4. Results 

First, test cases without rotation or magnetic diffusion were sim¬ 
ulated with field gradient p - 1. The Tayler instability is prop¬ 
erly traced for a number of growth timescales. More physical 
ingredients are built up each step in order to approach our sim¬ 
ple model with the conditions found in a radiative stellar interior. 

For the analysis of the simulations we need a formalism for 
measuring the amplitude of the various azimuthal m and vertical 
n modes, keeping in mind that n is the dimensionless vertical 
wavenumber; the n = 1 mode corresponds to a wavelength equal 
to the height of the computational box. We extract information 
from two fields: the change in the magnetic field bB = B - B(f = 
0), and the velocity field v. If we imagine the field to be frozen 
into the fluid, then bB reflects the displacement the integral of 
the velocity v. 

We can now decompose the component of these fields us¬ 
ing Fourier series in order to extract the amplitudes of the differ¬ 
ent azimuthal m modes. For instance with the velocity held 


f) = -Ao(?ir, z, f) + ^ ^R(A,„(tir,z,t)e""'^), 
;»=1 

where Aq and A^ are the Fourier coefficients given by 

1 r 

A()(ut, Z,t) - - (p, z, t) A(p 

and 

1 r 

A„,(ziT,z,t) - - I \(m,(j),z,t) cos(mtp)d)p. 

^ J-TT 

We then compute the amplitude of each m mode as 


(23) 


v_Amp,„it) = 


j-CTmax jL, 


vLvTT 


1/2 


(24) 


where tUnjax is set to cto/2, as the magnetic held strength in the 
region between zu - 0 to id-o/2 is approximately proportional 
to the distance from the axis to the power p, i.e. <x nrP. A 


similar procedure can be applied to the magnetic fields to derive 
the magnetic displacement amplitude, bB_Amp,„(f). 

Figure 6 shows the evolution of the amplitudes of the first 
three m modes as a function of time. The m - 1 mode is the only 
unstable mode, as predicted above in (5). When the amplitude 
of the unstable m - 1 mode is comparable to the wavelength, 
6B_Amp,„ ss the system goes into the non-linear regime. 
At this point energy can be transferred between modes and m + 1 
modes begin to grow. This is observed to happen around //ta ~ 
10 (where ta = I/wa) in figure 6. 

The displacement amplitude computed using equation (24) 
contains the contribution of all vertical wavelengths. For further 
analysis, it will be necessary to measure the amplitudes of each 
individual vertical n modes separately. A second Fourier decom¬ 
position, similar to the one applied in equation (23), is performed 
in the vertical direction. 


j /X—CXJ 

bv(nT, 0, z, f) = -C„,,o {rn, 0 + ^ ^R(C,„,„(nT, f) e'""), 

n=\ 


(25) 


where C„,_o and C,„^„ are the new Fourier coefficients given by: 

'-'L- nn 


= f 

Jo J-n 


and 




bv cos(ot0) d<p dz 


bv cos(ot0) cos(nz) di^ dz 


and are used to compute the displacement amplitude for each in¬ 
dividual pair of azimuthal m and vertical n modes, v_Amp^ Jf), 
similar to the procedure applied in equation 24. 


4.1. Case 1 - rotating diffusive star 

We first simulated the effects of just rotation and just diffusion 
on their own, for the constant-Alfven-frequency case, p = 1. 
Rotation is parallel to the magnetic axis. We compare our re¬ 
sults to those presented by Braithwaite (2006b) and show that 
we obtained similar results. 

^ Up to 8 or 9 orders of magnitude assuming a IG toroidal field in the 
Sun 
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Fig. 7. Amplitude 6B{t) of the m = 1, n = 1 mode as a function of 
time for various rotation rates. Blue lines correspond to rotation rates 
below the Alfven frequency, and red lines correspond to rotation rates 
above the Alfven frequency. 



Fig. 8. Growth rate of the Tayler instability for the ;« = 1 and n = 1 
mode as a function of the relative magnetic diffusivity rjlrjcnt measured 
from the velocity field, v_Amp^^{t). 


The rotation-only case has a threshold for the rotation rate 
equal to the Alfven frequency, Qy = wa,^- Rotation rates faster 
than that completely suppress the Tayler instability, presumably 
imposing some epicyclic motion of the perturbed fluid elements 
around an equilibrium point. Figure 7 shows the amplitude of 
the m- \,n - 1 mode in simulations with different fl||, ranging 
from Q||/wa,,^ = 0 - 10. It is clearly observed how the system 
is unstable for Q||/wa ,0 < 1 (blue lines), and is stabilized for 
rotation rates f 2 ||/wA,(S > 1 (red lines). 

In contrast, the diffusion-only case shows no threshold beyond 
which the Tayler instability is completely stabilized. 

Figure 8 shows the measured unstable growth rates in terms 
of iiri. It is observed that even at diffusivities as high as - 50; 
unstable growth is still present. 

For the simultaneous action of both rotation and magnetic 
diffusion a total of 56 simulations were conducted. Figure 9 
summarizes the measured growth rates normalized to the Alfven 
frequency, ct/wa.^. 

We consider a system to be stable if no unstable growth is 
observed for a hundred Alfven crossing timescales, IOOta, cor¬ 
responding to the white region in figure 9. We see that at very 
low magnetic diffusivities, left side of the figure, the system be- 





0.1 1.0 10.0 


Magnetic diffusivity B/Vcit 

Fig. 9. Contours of mean growth rate plotted as a function of nor¬ 
malized rotation rate tin = G||/tUA,,j and the relative magnetic diffusion 
n, = Jz/Jlcrit- The white region represents stability, for at least IOOta. 
The color contours correspond to increasing growth rates of the insta¬ 
bility, where hlack corresponds to small and blue to fast growth rates. 
The yellow dashed line corresponds to the relation = 2n, -l- 1 sepa¬ 
rating the stable and unstable regions. 

haves similarly to the non-diffusive case, the system is stable for 
rotation rates Qy > Now if we move from the stable region 
in a straight horizontal line towards the right side of the figure 
we go towards higher magnetic diffusive environments and pen¬ 
etrate a region of unstable growth. This is a very interesting 
result as one would prima facie expect that two stabilizing pro¬ 
cesses, such as rotation and magnetic diffusion, to collaborate in 
suppressing the instability, but instead the presence of magnetic 
diffusion disturbs the stabilization introduced by rotation. 

Presumably this happens because the Coriolis force pushes a 
fluid element onto an epicyclic equilibrium motion; magnetic 
diffusion constantly removes energy from the magnetic field 
such that there is no longer sufficient energy to return the fluid 
element to an equilibrium epicycle, and it will slowly spiral out¬ 
wards. 

The dashed yellow line in figure 9 is a relation between 
and n,j. If we assume that the boundary dividing the stable and 
unstable regimes can be separated by a linear relation we obtain 

«n|| ^ -I- 1. (26) 

It is not clear why the threshold should be here. The factor of 2 
may come from the 2 in the Coriolis force itself. 

Alternatively the stability condition can be expressed as 

(27) 

This is the stability condition for fast rotating, diffusive systems 
as a function of the perturbed wavenumber k. The growth rate of 
the instability is in general a function of the Alfven frequency, 
rotation rate, and magnetic diffusion rate, a = cr(cuA, fly, ? 7 ^^) 
for the non-stratified case. 

4.2. Case 2 - perpendicular rotator 

We simulated a rotator perpendicular to the magnetic field axis of 
symmetry. The magnetic field is oriented in the z-axis, B = e^, 
and the rotation axis is oriented along the tit((^ = 0 )-direction 
fi = ^±{ip - 0)enj. As in the previous section, the field gra¬ 
dient is p = 1. We expect rotation to have a different effect on 
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Fig. 10. Growth rate of the Tayler instability, normalized to the Alfven 
frequency, as a function of the perpendicular rotation rate, normalized 
to the predicted critical rotation rate for each n mode. The lines approx¬ 
imately on top of each other show the dependence of the damping on 
wavenumber is as predicted. 



Fig. 11. Growth rate cr of the Tayler instability as a function of the 
solid-body rotation rate Qy, in the p = 2 case. Normalized local growth 
rates are given for three different distances from the axis, cr(m). The 
dashed line corresponds to the predicted growth rate of the Tayler in¬ 
stability for slow rotators, cr = cja^- The dotted-dashed line is the 
predicted behavior of the instability for fast rotators, a = tu^^/2n||. 


different vertical wavelengths, as given by equation 9. Figure 
10 shows the growth rates of the first three vertical n modes 
(all m - 1 ) as a function of the rotation rate normalized to 
the critical rotation for each n mode. We can observe how the 
growth rate of the instability is damped with increasing rotation 
rate and furthermore how this damping affects different vertical 
modes. It is however important to remember that the horizon¬ 
tal axis of the plot is normalized to the critical rotation rate for 
each vertical mode. This means that given a rotation rate of e.g. 
Qj^ = 4Qx,ciit(« = 1 ) = '2-0(jJA,<p, the n - \ mode is damped by 
~ 90%, n - 2 mode by ~ 45%, and n - mode by ~ 25%. In 
agreement with our prediction that perpendicular rotation damps 
different vertical modes with a different strength. 

Even if we did not see complete stability in our simulation, 
figure 10 suggests that there is stabilization for each vertical 
mode at a rotation rates of Qj^ ^ 5Qj.,crit(n). Apparently our 
stability condition (9) is missing a factor of order unity. This is 
not surprising, given that the damping effect should depend on 
the precise arrangement of the displacements within the mode. 


4.3. Case 3 - quadratic radial field gradient 


We saw above that in the case where p - I, rotation can com¬ 
pletely stabilize the magnetic field, at least where there is no 
diffusion. We now investigate the case with a radial field gra¬ 
dient of p = 2, which is predicted to behave differently. The 
Alfven frequency grows proportional to distance from the axis 
a>A,<t>{'nT) oc nr. For this reason we expect that the instability 
grows faster at larger radii than close to the axis. We introduce 
a new quantity, bZ?_Amp„, „(zit, f), corresponding to the displace¬ 
ment amplitude, for a given pair of azimuthal and vertical modes, 
as a function of time and distance to the axis nr; 


„(nr, t) 


/ C*,_„(nT, f)Cm,„(trr, t) 
\ nr“L.n 



(28) 


This new quantity allows us to study the growth rate of the insta¬ 
bility at different distances from the axis. 

Local growth rates, cr(nr), are computed from the local am¬ 
plitude displacements, bZ?_Amp,„,j(nT, f)- This quantity is then 
divided by the local Alfven frequency, oja,( pinr), and the normal¬ 
ized local growth rate is obtained, cr(nr)/niA, 0 (tir). 


Figure 11 shows the behavior of these normalized local 
growth rates as a function of rotation rate. The first main feature 
of the figure is that even for the very fast rotators, instability is 
still present. This is in agreement with the prediction that insta¬ 
bility is always present in fast rotating stars containing a purely 
toroidal field with a radial field gradient p > 3/2, obtained from 
the stability conditions ( 8 ). 

There are two behaviors observed in this plot separated by a 
critical rotation rate at Qy wa, 0 . In the slow rotation regime 
where Q||/wa ,0 < 1 , the growth rate cr cja,^, whereas in 
the fast rotation regime where Qy/cuA.^ ^ 1 , we have instead 
cr aj ^/2Qy. Which is exactly what Pitts & Tayler (1985) and 
Spruit (1999) predicted a couple of decades ago."^ 

4.4. Case 4 - mixed poloidal-toroidal field 

This section is intended to shed some light in the stability of 
mixed poloidal-toroidal magnetic field configurations in radia¬ 
tive stellar interiors. We focus our analysis on the behavior 
of the instability known to be present in purely toroidal fields, 
but now in mixed poloidal-toroidal field configurations. As dis¬ 
cussed in section 2.3 a poloidal field could damp, and if strong 
enough probably be able to suppress, the instability of a toroidal 
field. Tayler’s (1973) stability conditions however are not valid 
for mixed field configurations. 

We note that the results presented in this section should be 
reviewed carefully because of the appearance of an unknown, 
spurious mode in the simulations growing always at the smallest 
scale in the vertical direction^. Such a mode was absent in 
purely poloidal and purely toroidal field simulations. We work 
under the hypothesis that any well resolved unstable mode 
growing faster than this unknown perturbation is real, and 
relevant for our analysis. The nature of the unknown mode 

Except that they neglected the factor of 2. At a latitude of 30" this 
factor will go away anyway since the component of fJ perpendicular 
to the spherical surfaces in which the fluid is moving around is half as 
much as at the poles. 

^ Despite the many vertical resolutions tested: 8, 16, 32, 64 and 96 
these perturbations, with vertical wavelength corresponding to that of 
the Nyquist frequency, was present in the m = 1 mode with a growth 
rate » (jJa/2. 
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field (5) and it is being perturbed in the mixed configuration. 
Growth of the m - 0 mode begins to appear with increasing 
poloidal field strength, reaching a maximum near the critical 
poloidal-to-toroidal field ratio, P„ ^ \. 

The presence of this unstable mode, not predicted analyti¬ 
cally, might be important for the hydromagnetic dynamo theory, 
given that this instability triggers axisymmetric vertical motions 
which it is often claimed are necessary for closure of the dynamo 
loop. We will investigate the importance of this mode in a future 
paper. 


Fig. 12. Growth rate of the first two azimuthal m modes as a function 
of the pitch-numher. The solid black line corresponds to the m = 0 
mode, stable at very high and low pitch-numbers but with a peak in 
unstable growth near the critical pitch-number, P„ ^ 1. The red dashed 
line corresponds to the m = I mode, unstable for very high pitch- 
numbers and stabilized near the critical pitch-number, k 1. 


will be studied in detail in a forthcoming paper and it must be 
established whether it is astrophysically relevant. 


We first introduce some important definitions to avoid mis¬ 
understandings; the azimuthal Alfven frequency, = VA.^/trr, 
is calculated with the toroidal Alfven velocity. In the initial con¬ 
ditions, the radial field gradient of the toroidal field is chosen 
for simplicity to grow linearly with the distance from the axis, 
p - I, therefore the Alfven frequency is constant throughout the 
volume of interest. The vertical Alfven frequency cja.z is calcu¬ 
lated with the vertical Alfven velocity, Va,z and is given by 


VA,; ^ B, 
L: Lj 


(29) 


The poloidal field strength is constant everywhere, - Bq, thus 
the vertical Alfven velocity and the vertical Alfven frequency 
are constant. 


For this field configuration field lines follow spirals as shown 
in figure 5b. The number of spirals a field line completes in the 
box depends on the relative strengths of the toroidal and poloidal 
fields. The more complete spirals a field line draws, the stronger 
the toroidal field is. We call this number of spirals in the box the 
pitch-number, P„, which is given by 


Wa,0 
2nu)A,z' 


(30) 


re-writing the instability condition (14) in terms of the pitch- 
number and solving for the « = 1 mode, we obtain that instability 
is expected to be present for > 1. We note that is only 
independent of nr in the special case p - \, which is what we 
use here. 

Figure 12 shows the behavior of the first two azimuthal 
modes, m - 0 and m - 1, as a function of pitch-number. 
The m - 1 mode has a growth rate on the order of the Alfven 
frequency for a predominantly toroidal field configuration, i.e. 
large pitch-number, and drops suddenly when the pitch-number 
reaches a critical value. This drop of the growth rate agrees with 
the stability conditions (14). We note that no complete suppres¬ 
sion of the instability is observed, even for very small pitch- 
numbers. However, these small growth rates are the effect of 
non-linear interaction with the spurious modes. 

The second interesting mode shown in figure 12 is the m - 0 
mode. This mode was marginally stable for the purely toroidal 


5. Discussion and summary 

Solid-body rotation, magnetic diffusion and field configuration 
are expected to have a strong infiuence on the presence and evo¬ 
lution of hydro-magnetic instabilities in toroidal fields in radia¬ 
tive stellar interiors. We present a set of 3D MHD simulations 
of purely toroidal and of mixed poloidal-toroidal magnetic fields 
in non-convective stellar interiors. A combination of various ra¬ 
dial field strength gradients {p — dlnB/dlnnr), solid body rota¬ 
tion rates (Q), magnetic diffusivities ( 77 ), and poloidal-to-toroidal 
field strengths (F„) is scoped. 

We observe the presence of the kink type m - 1 Tayler 
instability in purely toroidal fields. We measure its growth 
rate, and compare our results with the various analytical results 
for non-rotating and non-diffusive (Tayler 1973), fast rotating 
(Acheson 1978; Pitts & Tayler 1985), highly diffusive (Spruit 
1999), and rotating-diffusive (Pitts & Tayler 1985; Spruit 1999) 
radiative stellar interiors. We also study the development of 
instabilities in a mixed poloidal-toroidal field configuration. 

Magnetic diffusion is not able to suppress the occurrence of 
the instability even at high diffusivity, in accordance with the 
results of Braithwaite (2006a). 

Diffusion and rotation parallel to magnetic axis: If the 
magnetic and rotation axes are parallel, fast rotation (G|| > wa,^) 
is able to stabilize the configuration (see figure 7 ), where aiA.fS 
is the azimuthal Alfven frequency. This is in the case where the 
field gradient p - However as soon as magnetic diffusion is 
included, the stability conditions now depend on a combination 
of the rotation rate, the magnetic field strength and the magnetic 
diffusion rate given by equation 27, Which suggests that very fast 
rotation can stabilize the system only if the magnetic diffusion 
rate does not subtract energy fast enough such that Coriolis force 
is unable maintain the fluid elements in stable epicyclic motions. 

Perpendicular rotation: If the magnetic field axis and the 
rotation axis are perpendicular, the damping strength introduced 
by rotation is different for different vertical wavelengths. The 
damping effect is stronger for small vertical wavenumbers k 
compared to large wavenumbers because of the larger defor¬ 
mation occurring perpendicular to the rotation axis at smaller 
wavenumbers. It is found that the instability is suppressed if 
Qj_ > 5wA,(StJTok, where mo is the length scale of the motions in 
the horizontal direction. 

Radial field gradient: To analyze the effects of the radial 
field gradient and rotation in the instability we investigated two 
cases, the linear radial field gradient, p - 1 and the quadratic ra¬ 
dial field gradient, p - 2. As described above, in the p - 1 case 
(without diffusion) stability is reached for fast rotation (meaning 
faster than the azimuthal Alfven frequency). In this p - 2 case, 
very fast rotation is never able to stabilize the configuration. Lin¬ 
ear growth rates behave broadly as predicted by Pitts & Tayler 
(1985), where for slow rotators the growth rate is roughly equal 
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to the Alfven frequency, whilst for fast rotators, the growth rate 
is damped by a factor of wa,0 /2Q||. It is probably safe to as¬ 
sume that this is also the case for all values of p above 3/2, as 
originally predicted. We have therefore provided a check on the 
theory, useful in light of the fact that the originators of the ana¬ 
lytic results were, we think it is fair to say, less than 100% certain 
they were correct. 

Mixed poloidal-toroidal field: Finally we simulated a 
mixed poloidal-toroidal field, and in particular the stabilizing 
influence of a uniform poloidal field. In the p - \ case, where 
in the absence of a poloidal held, only the m - \ mode is 
unstable, we And that the m - I mode is indeed suppressed by 
the poloidal held, and that the stability criterion is where it is 
predicted to be. Interestingly, we also And that the m - 0 mode 
is unstable in the vicinity of this threshold. 

The presence of the unstable axisymmetric m - 0 mode 
in the mixed poloidal-toroidal field configuration may be quite 
important. In the case where the field gradient is such that 
it is marginally stable in the purely toroidal case, it becomes 
dominant in the vicinity of the threshold of stabilization by the 
poloidal field, which of course is precisely where the dynamo is 
likely to be. It opens the possibility for axisymmetric turbulent 
radial motions which affect the mixing and the regeneration of 
the poloidal field. It is often claimed by various members of the 
MHD community that a poloidal field cannot be renewed with¬ 
out the m = 0 mode; whether this claim is correct is, if the m = 0 
mode is present and perhaps even dominant, no longer relevant. 

Indeed, as far as the field gradient p is concerned, in the 
fast-rotating (which is expected to occur much more frequently 
in nature than the slowly rotating) adiabatic case we can sum¬ 
maries behavior in terms of the threshold; where p < 1 all 
modes are stable; where 1 < p < 3/2 just the m = 0 mode 
is unstable; where /? > 3/2 the m = 1 mode is also unstable, 
with thresholds at higher values of p introducing instability at 
progressively higher values of m. However, these thresholds are 
probably of little practical importance. Whilst it is often said 
that the Tayler instability is dominated by the m = 1 mode be¬ 
cause it is that mode which is in some sense the first to become 
unstable, this is misleading: only in the rather special circum¬ 
stance of -1/2 < /? < 1 and Qy < wa is the m - 1 mode the 
only one present. Firstly, the second of those two conditions is 
rarely fulfilled in nature. Secondly, in the context of the fluc¬ 
tuating field envisaged in the Tayler-Spruit dynamo model, the 
field gradient p will have values all the way from p = -oo to 
p = - 1 - 00 . It seems perfectly plausible that the dynamo could be 
dominated by the m - 0 mode, or indeed by m = 20. Roughly 
speaking, the instability is present where p is positive and absent 
where it is negative: the dynamo can therefore be expected to be 
intermittent in some sense. 
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